Aging and fluctuation-dissipation ratio for the diluted Ising Model 
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We consider the out-of-equilibrium, purely relaxational dynamics of a weakly diluted Ising model 
in the aging regime at criticality. We derive at first order in a y/e expansion the two-time response 
and correlation functions for vanishing momenta. The long-time limit of the critical fluctuation- 
dissipation ratio is computed at the same order in perturbation theory. 
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According to universality hypothesis, critical phenom- 
ena can be described in terms of quantities that do not 
depend on the microscopic details of the systems, but 
only on global properties such as symmetries, dimen- 
sionality etc. A question of theoretical and experimental 
interest is whether and how this critical behavior is al- 
tered by introducing in the systems a small amount of 
uncorrelated impurities leading to models with quenched 
disorder. 

The static critical behavior of these systems is well un- 
derstood thank to the Harris criterion [jj]. It states that 
the addition of impurities to a system which undergoes a 
second-order phase transition does not change the critical 
behavior if the specific-heat critical exponent a p of the 
pure system is negative. If a p is positive, the transition 
is altered. 

For the very important class of the three-dimensional 
0(M)-vector models it is known that a p < for M > 
2 [jj , and the critical behavior is unchanged in presence of 
weak quenched disorder. Instead, the specific-heat expo- 
nent of the three-dimensional Ising model is positive [^) , 
thus the existence of a new Random Ising Model (RIM) 
universality class is expected, as confirmed by Renormal- 
ization Group (RG) analyses, Monte Carlo simulations 
(MCs), and experimental works (see Refs. for a com- 
prehensive review on the subject, and for an updated list 
of references). 

The purely relaxational equilibrium dynamics (Model 
A of Ref . Q ) of this new universality class is under inten- 
sive investigation [^]-^| . The dynamic critical exponent z 
differs from the mean-field value already in one-loop ap- 
proximation ||, at variance with the pure model. This 
exponent is known at three-loop level in an yfe [Q] and in 
fixed (d = 2, 3) dimension [|| expansion, and has a value 
in good agreement with several MCs ||. 

The out-of-equilibrium dynamics is instead less stud- 
ied. The initial slip exponent 9 of the response function 
was determined at two-loop order [ [To[ and the response 
function only at one-loop, both for conservative and non- 
conservative dynamics Jll| . 

In this work we focus on a different regime of out- 
of-equilibrium dynamics: the aging one. The relaxation 
of the system from an out-of-equilibrium initial state is 



characterized by two different regimes: a transient one 
with off-equilibrium evolution , for t < and a sta- 
tionary one with equilibrium evolution of fluctuations for 
t > £r, where t R is the relaxation time. In the former 
a dependence of the behavior of the system on initial 
condition is expected, while in the latter homogeneity of 
time and time reversal symmetry (at least in the absence 
of external fields) are recovered. Consider the system in 
a disordered state for the initial time t = 0, and quench 
it to a given temperature T >T Cl where T c is the critical 
temperature. Calling 4> x (t) the order parameter of the 
model, its response to an external field h applied at a 
time s > and in x = is given by the response func- 
tion R x (t, s) = 8((j) x (t)) /5h(s), where (• • •) stands for the 
mean over stochastic dynamics. The two-time correlation 
function C x (i, s) = (0 x (£)0o(s)) is useful to describe the 
dynamics of order parameter fluctuations. 

When the system does not reach the equilibrium (i. e. 
ifl, = oo) all the previous functions will depend both on 
s (the "age" of the system) and t. This behavior is usu- 
ally referred to as aging and was first noted in spin glass 
systems ]l2]] . To characterize the distance from equilib- 
rium of an aging system, evolving at a fixed temperature 
T, the fluctuation-dissipation ratio (FDR) is usually in- 
troduced [pl: 



X x (t,s) 



TR x {t,s) 



(1) 



d s C x (t, s) 

When t,s 3> tn the fluctuation-dissipation theorem holds 
and thus A x (i, s) = 1. 

Recently Jf3|-pcH attention has been paid to the FDR, 
for nonequilibrium and nonglassy systems quenched at 
their critical temperature T c from an initial disordered 
state | pl[ . The scaling form for i? x= o was rigorously 
established using conformal invariance p0| , p2[ |. Within 
the field-theoretical approach to critical dynamics, cal- 
culations are simpler if done in momentum space, thus 
we are interested in momentum-dependent response and 
correlation functions. From RG arguments it is expected 
that they scale, for q = 0, as ||,|w|,|l|,^9| 

T c i? q=0 (t, s) = A R (t - S ) a {t/ S ) e F R { S /t) , (2) 
C q=0 (t, s) = A c s{t - s) a (t/s) 9 F c ( s /t) , (3) 
dsC^oit, s) = A 9c (t - s ) a (t/ s ) e F dc ( s /t) , (4) 
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where a = (2 — r\ — z)j z. The functions Fc{v), Fgc(v) 
and Fr(v) are universal provided one fixes the nonuni- 
versal normalization constant Ar, Ac, and Aqc to have 
Fi(0) — 0. Obviously Ac, Aqc, Fc(v), and Fgc(v) are 
not independent, in fact it holds 



A dc = A c (l - t 
F dc (v) = F c (v) + 



1 



F' c {v) 



1-v 



F c (v) 



(5) 
(6) 



It has been argued that the limit X^ Q = 
limbec , limj—xx, X x=0 (t, s) is a novel universal quantity 
of noncquilibrium critical dynamics [^5|jr^ 1 and may be 
computed as an amplitude ratio. Also the analog in q 
space |lq| 



*q=o(*, S) 



A, 



A j 



OBq=o(M) _ 

d s C q=0 (M) A dC A c (l-e) 



(7) 



is universal. In Ref. Jig ] an heuristic argument is given 
to show that X™ =Q = X£L Q . 

In a recent work |l9| we evaluated X^ at critical- 
ity for 0(M)-vector models at the second order in the 
e-expansion, finding values in very good agreement with 
two- and three- dimensional numerical simulations for the 
Ising model |l5 ,[l7j . We also confirmed the validity of the 
scaling laws (p) and (0) at the same order in perturba- 
tion theory. A different model of relaxation (Model C of 
Ref. [§) has also been studied @. 

The extension of this kind of investigation to disor- 
dered systems is very interesting because, besides giving 
a check of the expected scaling laws, it predicts a new 
universal dynamical quantity (the long time limit of the 
FDR) which could be measured in MCs and could be 
used to identify a universality class, as in the case of 
other universal quantities. 

The time evolution of an A-component field <£>(x, t) 
under a purely dissipative relaxation dynamics (Model 
A of Ref. 0|) is described by the stochastic Langevin 
equation 



d t <p(x, t) 



(8) 



where f2 is the kinetic coefficient, £(x, t) a zero-mean 
stochastic Gaussian noise with correlations 

fe(x,^ i (x',t')) = 2^5(x-x')^-0^- (9) 
and H.^[ip] the static Landau-Ginzburg Hamiltonian || 

H» = J ddx \{dvf + ^(ro + ^(x))^ 2 + 

(10) 

Here V'( x ) is a spatially uncorrelated random field with 
Gaussian distribution 



P(V) = 



: exp 



Aw 



(11) 



Dynamical correlation functions, generated 
by Langevin equation (||) and averaged over the noise 
£, can be obtained by the field-theoretical action pa] 



S^,0\ = dtd d 



.dip n .dH^[ip] 



<piltp 



(12) 



where y>(x, t) is the response field. 

The effect of a macroscopic initial condition <ySo( x ) = 
tp(x,t — 0) may be taken into account by averaging over 
the initial configuration with a weight e^ H °^ Vo ', where 



x Y^o(x) 



a(x)) 2 



(13) 



This specifies an initial state a(x) with short-range cor- 
relations proportional to t^ 1 . 

In this way all response and correlation functions 
may be obtained as averages on the functional weight 
exp{— (S^Yp, 0\ + Ho[ipo])}. In the analysis of static crit- 
ical behavior, the average over the quenched disorder ip 
is usually performed by means of the replica trick j^] . If 
instead we are interested in dynamic processes it is sim- 
pler to perform directly the average at the beginning of 
the calculation Eq] 



[dV]PW)exp(-^M) = exp(-S[<p,fl) (14) 



with the -0-independent action 



s[ Y - . ;•] = /<i''./ <J j dt<f[d t (p + n(r -A)tp-riLp} 



-irj dt ^~ 2 



d*W h ( 15 ) 



where v$ oc w. 

The perturbative expansion is performed in terms of 
the two fourth-order couplings go an d vq and using the 
propagators of the free theory with an initial condi- 
tion at t = 0, (</5,(q, s)<pj(— q,i))o = <%i?°(£, s) and 
((Pi(q,s)(pj(-q,t)) = SijC°(t,s) |||. It has also been 
shown that Tq -1 is irrelevant (in RG sense) for large times 
behavior p3"[ . From the technical point, the breaking of 
homogeneity in time gives rise to some problems in the 
renormalization procedure in terms of one-particle irre- 
ducible correlation functions (see Ref. |2jJ and references 
therein) so all the computations are done in terms of con- 
nected functions. 

To compute the response function at one-loop level, we 
have to evaluate the two Feynman diagrams depicted in 
Fig. [l]. In terms of them we may write 

Rq(t,s) =R°(t,s)- -go(a) + v (b) + 0(gl,v$, g v Q ) , 

(16) 
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FIG. 1. Feynman diagrams contributing to the one- loop 
response function. Response functions are drawn as 
wavy-normal lines, whereas correlators are normal lines. A 
wavy line is attached to the response field and a normal one 
to the order parameter. The dotted line is a non-local v-like 
vertex. 



where we are considering the case N = 1 (RIM univer- 
sality class), and we set 17 = 1 to lighten the notation. 

In the following we report the expressions of Feynman 
diagrams at criticality (tq — in dimensional regulariza- 
tion) for vanishing external momentum, since we are only 
interested in that limit, and since expressions for nonzero 
q are long and not very illuminating. 

The diagram (a) in Fig. [l] contributes also the the re- 
sponse function of nondisordered models, and it has been 
computed in Ref. obtaining (for t > s): 



(a) = -N d ]\og- + 0(e). 
4 s 



(17) 



where N d = 2/[r(d/2)(47r) d/2 ]. For diagram (b) we find 



(h) j d/'d/" l-^Ro(t,t')R p (t',t")R (t",s) 



1 1 1 

( 47r )d/2 i _ d / 2 2 - d/2 



(t-s) 



2-d/2 



(18) 



Inserting the expression for (a) and (b) in Eq. (JO]) and 
expanding (b) at first order in e, one obtains 



R B (t,s) = 1 



- 1 ^ 1 
.9o o m_ 



2 



log(t 



■IE 



-0(e ,eg ,ev ,g ,v ,g v ) , 



(19) 



where go — Ndgo and vq = NdVQ. The dimensional pole 
in this expression can be canceled by a multiplicative 
renormalization of both the fields tp and (p. 

The critical response function is then obtained setting 
the renormalized couplings at their fixed point values. 
We remind that the stable fixed point of the RIM is of 
order ^fe and not e (see, e. g. Q), due to the degeneracy 
of the one-loop f3 functions. The nontrivial fixed point 
values at the first non-vanishing order (i. e. two loops) 
are 



,9o = .9 



+ 0(e). (20) 



Finally we get 
R(t,s) = 1 + 



In- 



ln(t - s) - j E 



0(e) 



(21) 



that is fully compatible with the expected scaling 
form (||) with the well known exponents [ ?],[ll) 

1 Rk s „ 1 /6e 



0(e) 



2 V 53 w ' 2 V 53 

and the new results Fr(x) = 1 + 0(e) and 

R = 2V 53 7B+ W ■ 



+ 0(e), (22) 



(23) 



There are five diagrams contributing to the correla- 
tion function. Four of them are obtained by the ones of 
Fig. ^ changing one of the two external response propa- 
gators with a correlation line (see Ref. jl9| for a detailed 
explanation of this correspondence). We call these four 
diagrams (a±), (02), (bi), and (62)- The sum (ax) + (02) 
was computed in |18) leading to 



(ax) + (a 2 ) 



N d ( t 
■s I log ■ 



2 V - s 
The sum (bx) + (62) is instead 

N d T(d/2) 



0(e). 



(24) 



(bx) + (b 2 )= 



(l-d/2)(2-d/2)(3-d/2) 

[t3-d/2 + s 3-d/2_ (t _ s) 3-d/2] 



(25) 



The octopus diagram in Fig || does not have a corre- 
sponding one contributing to the response function. It 
has the value 



(c) 



N d T(d/2) 



(l-d/2)(2-d/2)(3-d/2) 

(t 8?-W + (t + , d/2 ,_ d/2 



(26) 



Collecting together these contributions and expanding in 
powers of e we find, at 0(e 2 , eg , ev , go,v%, g v ), 

C B (t, s) = 2s- ^((ax) + (02)) + v ((bx) + (b 2 ) + (c)) 



= 2s 



.9o (, t 

T s log 7 



2s 

lES 

e 



(27) 



(t - s) log(* - s) (t + s) log(< + s) 




FIG. 2. Feynman diagram contributing to the one-loop cor- 
relation function that does not have analog in the response. 
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Also this expression is renormalized by canceling the di- 
mensional pole by means of multiplicative renormaliza- 
tions of parameters and fields. At the nontrivial fixed 
point for couplings we have 



C(t,s) 




log - + 3 - 7£ - log(t - s) 
s 



t + s 



0(e) 



(28) 



which is again compatible with the expected scaling form 
(cf. Eq. (||)) with the same exponents of Eq. p2]), the 
non universal amplitude 



Ac 



= 1 



(2 - y E ) + 0(e) , 
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and the universal regular scaling function 



F c (x) 



1 



1 



log- 



1 



(29) 

0(e). 
(30) 



Note that at variance with the pure model [[18 19 1, the 
function Fc(x) receives a contribution already at one- 
loop order which should be observable in MCs. 
Using Eq. @) the FDR is 



X*=« = \-l\rT- o +0(e) 



(31) 



that, for e = 1 leads to X£L ~ 0.416, and - 0.381 for 
e = 2. It would be interesting to see if this one- loop 
result is in as good agreement with MCs as in the case 
of the pure model (cf. Refs. Jl^,[ll|). To this order it 
is not even clear whether randomness really changes in 
a sensible way the limit of the FDR or not. Two- loop 
computations and MCs could clarify this point. 

For completeness we report also the FDR for finite 
times: 



1, t-s 
— log 

2 5 t + s 



0(e). (32) 
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